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We propose a method to probe the equation of state of the early universe and its evolution, 
using the stochastic gravitational wave background from inflation. A small deviation from 
purely radiation dominated universe (w = 1/3) would be clearly imprinted on the gravitational 
wave spectrum ^lcw(f) due to the nearly scale invariant nature of inflationary generated waves. 

7— I . 

0^ . Gravitational waves generated in the early universe are one of the most interesting targets in 

O ■ 

l/~s ■ observational cosmology. Though direct detection of waves is not a simple task due to their ex- 
■ tremely weak interaction, they could serve as an invaluable fossil of the early universe that can be 
hardly probed by other observational means. 1 ' 2 ) The most realistic generation mechanism of the 

O ! . . . . . .ON 

q-i stochastic gravitational wave background in the early universe is quantum effects during inflation. ' 

While the fact that we live in a large and nearly homogeneous universe alone is sufficient to prove 
the existence of an inflationary epoch in the early universe, ' recent observations of the WMAP 



satellite 5 ^ have provided a more sophisticated confirmation of inflation. This comes from that 
super-horizon-scale perturbations that were generated in the very early universe present a negative 
correlation between temperature anisotropy and E-mode polarization (TE). 6 ) Note that observa- 
tion of nearly scale-invariant spectrum of curvature perturbation alone is insufficient as a proof of 
inflation because other causal mechanisms can also exist to generate a similar spectrum. 7 ) Such 
mechanisms, however, predict a positive correlation in TE and is now discriminated by WMAP 
observations. 6 ) 

Long wavelength components of the stochastic gravitational wave background generated during 
inflation can in principle be detected by the B-mode polarization of cosmic microwave background 
radiation (CMB), 8 ) and shorter- wave components with frequency ~ Hz may be observed by a space 
laser interferometer. As for the former, although WMAP has obtained only an upper bound, 9 ) a 
number of observational projects are planned, and those include Planck 10 ) and BICEP 11 ) suggesting 
a strong possibility for detection of the stochastic background wave from inflation by CMB. On the 
other hand, most theoretical studies of the gravitational wave background of shorter wavelength 
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have emphasized on the detection itself so far. However a discussion with an anticipation for future 
diverse possibilities of gravitational wave cosmology, beyond the first detection, is worthwhile. 

In the present paper, we argue that the amplitude of high-frequency stochastic gravitational 
wave background not only carries information on the inflationary regime, during which they are 
generated, but also serves as a probe of the equation of state in the early universe through which the 
waves propagate towards us. This double role played by primordial gravitational wave background 
in particle cosmology is similar to that played by high-redshift quasars in observational cosmology, 
which not only reflects the properties of high-redshift universe, at their location, but also carry 
line-of-sight information in their absorption spectra. In this sense, just as high-redshift quasars are 
regarded as lighthouses of the distant universe, we may regard the inflation-produced gravitational 
wave background as a lighthouse that can shed light on the early universe physics. 

As is well-known, equation of state of the early universe, just after t = 1 sec and T = 1 MeV, 
has been accurately probed by big-bang nucleosynthesis (BBN) calculations and a comparison of 
their yield with observations, 12 ** though we do not have any sensible means to probe it prior to 
this epoch. This is why we only impose a mild constraint that the universe should be dominated 
by radiation before BBN when we work on particle-physics model building of the early universe 
involving lately decaying particles. Thus, if we could observe the spectrum of gravitational waves 
from a source with a known property in the early universe, we could obtain important information 
and constrain high energy physics with which we describe the evolution of the universe. 

We start with a brief review of properties of a gravitational wave generated quantum mechanically 
during inflation with the purpose of demonstrating our notation. We introduce tensor perturba- 
tions, hij, around a spatially flat Robertson- Walker metric as 

ds 2 = -dt 2 + a 2 {t) (Sij + 2hij) da? do?, (1) 

with a being the scale factor. Decomposing the tensor metric perturbation to Fourier modes as 

ha = V^G £ J j^Lfii^etj, (2) 



A=+,x ' 



we find that the two independent degrees of freedom ip behave as two massless minimally coupled 
scalar fields, where represents polarization tensor with e^-e'^' = 5 AA ' for A, A' = +, x. 
Applying quantum field theory of a massless minimally coupled field in de Sitter spacetime, we find 
that the Fourier modes are characterized by the following vacuum correlation, 

(^{t)^;{t)) = ^{k-k')8 AA \ (3) 

so that the amplitude per logarithmic frequency interval is given by 

*l</> = *<M«</)Hx^(^§(^)', (4) 
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where Mp\ = G~ 1 / 2 denotes the Planck scale. Here H((j>) denotes the Hubble parameter during 
inflation to be evaluated when the mode with comoving frequency / left the Hubble radius during 
inflation. It is expressed as 

in terms of the inflaton potential V[4>], where <fi(f) denotes the value of the inflaton field when the 
tensor mode / crosses the horizon during the inflationary epoch. 

In this paper, we use the present frequency / to characterize the comoving wavelength of gravi- 
tational waves, and the scale factor a for the cosmic time t. For the spectrum of gravitational wave 
background, it is convenient to use a dimensionless quantity &Gw(f, a) defined by 

n (f x PG\v(f,a) s _ h 2 {f,a) /2vr/a \ 2 . . 

^ (/ ' a) = ^^T' pGw{f ' a) = ^MG\—) ' (6) 

where pcw{fi a ) is the energy density of gravitational waves per unit logarithmic frequency interval 
around /, ao is the current value of the scale factor, and p C r{a) is the critical density of the universe. 
In the early universe we can identify the critical density p C r{a) to the total density of the universe 
Ptot(d), as the spatial curvature is negligible. 
In a power-law background a(t) oc t p with p < 1, each Fourier mode behaves as 

i-3p ( V k \ 

h(f, a) oc a(t) J 3 P -i ,,,„,,, , k = 2vr/a , (7) 

2(i-p) V 1 — P a(t)H (t) J 

with Bessel function J n (x). Such that the amplitude of gravitational wave takes a constant value, 
h(f,a) = hp{f), until a mode reenters the Hubble radius H~ l at a = 2itfao/H = aj n (/), where 
the normalized energy density reads 

Evolution of the relative energy density of the tensor modes within the Hubble horizon depends 
on the equation of state of the universe as we see below. From the definition of VLgw, we have 

dlnn G w{f, a) = din pcwjf, a) _ din p to t{f, a) 
din a din a din a 

From eqs. (6) and (7), the tensor modes well within the horizon behave simply as radiation, 

din pGwjf , a) 
din a 

and evolution of the total energy density is given by 

d\np tot (f,a) , n] , . 

=-31 + «/o , (11) 

a ma 

where w is defined by w = P/ptot with pressure of the universe, P. From eqs. (9-11) we have 

dlnn G w(f, a) „ ( v 1 no , 

— = 3w(a) - 1, 12 

a ma 



(9) 



= -4, (10) 
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or 

raa 



raa 

\a.Q GW (f,ao) = [3w(a) - l]d\na + lnQ GW (f,a in (f)). (13) 

For our analysis it is convenient to use the following expression that relates the present energy 
density Q G \y(f,ao) at different frequencies f± and f 2 , 

ra-inih) V\6(fr>)] 

\nn GW (f 2 ,a )-\nn GW (f 1 ,a )= [3w(a) - 1] din a + In -g ) J2 . (14) 

The first term represents the effect caused by deviation of the equation of state from w = 1/3 
and the latter one shows the difference of the intrinsic amplitudes £l G w(f, o-in(f)) generated at the 
inflationary epoch. Using the slow-roll equation of motion, 

3H4> = -V'[4>], (15) 

we can write down the latter term as 

lr V[0(/i)] " Svr \V) m l/J' Ubj 

with = Here the coefficient M pl /8ir (V /V) 2 is evaluated around 4> ~ 0(/i) and 

would be much smaller than unity, which is one of the slow-roll parameters. 

As is well known, under the assumption that fluctuations are governed by a single inflaton field 
cf), we can obtain the magnitude and the first derivative of the inflaton potential V[cf>] from the 
quadrupole anisotropy of the cosmic microwave background radiation (CMB). Denoting the scalar 
and the tensor contributions by S and T, respectively, we find 

_ 2.2{V/M* l ) / V \ 

Sl=2 ~ ( Mpi v/vr 1=2 ~ kkJ' { ) 

where V or V' is evaluated at 4>(f) with mode / corresponding to the present horizon 2irf = Hq, and 
the values of the numerical coefficients are those in the case of the Einstein de Sitter universe. 13 ' 14 ) 
As for the spectral indices ns and ht for these two modes, there are deviations from the results 
(ns = 1 and ut = 0) for purely de-Sitter inflation due to variation of the potential V [</>(/)] along 
with the modes /, 13 ' 14 - ) 

M 2 Pl (V'\ 2 M 2 Pl (V'\' _M 2 Pl fV^ 2 



With eqs. (17) and (18) the coefficient M pl /8ir (V' /V) 2 is expressed in a well-known form 

Mpi ( v '\ 2 r - T i=2 no , 

77 = n T = ~^ r =77-- ( 19 ) 



8vr \V J 7.0' S l=2 

However, the factor 7.0 in the above expression should be replaced by 5.0 for realistic values of the 
cosmological parameters in the concordance model. 15 ) WMAP has obtained an upper bound on the 
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coefficient in eq. (16), as |Mj, z /87r (V'/V) 2 | ~ r/5 < O.l. 9 ) Thus even in the frequency difference 
of an order of magnitude, ln(/2//i) ~ 2, the second term of eq. (14) would be very small. If we 
could measure the difference of lnQcw(f2, o-o) — ^ n ^Gw(fii ^o), for example to be ~ 0.5, it would 
mainly be due to the first term. Namely, we could detect a curious behavior of equation of state 
w 7^ 1/3 in the early universe. 

As argued above, there are essentially two frequency bands which may be regarded as a serious 
target for observations. One is the ultra-low frequency band with / ~ H , which can be detected 
indirectly through B-mode polarization of CMB. There are some CMB projects aiming to measure 
the tensor fluctuations as a primary target. As explained earlier, BICEP, which is expected to start 
observations in 2004, has a sensitivity down to r ~ O.l. 11 ) CMBPOL, a space mission planned to 
be launched around 2012, has sensitivity down to r ~ 0.0001 which is close to the fundamental 
limitation for the detection of r using CMB polarization. 16 ) Planck, which is scheduled to be 
launched in 2007, 10 ) will determine the index ns with error less than An,s < 0.005. With these 
future CMB projects, we will observationally obtain information related to the potential V[<j>\ at a 
large scale, / ~ Hq, and can use it to estimate V[<p] up to the second derivative using (18). Thus, 
although this frequency band is irrelevant to the evolution of the cosmic equation of state, the 
information obtained through these observations are invaluable that allow us to extrapolate the 
spectrum £lGw(f,ain(f)) up to much higher frequencies 14 ) which can be used as a calibrator in 
order to investigate an anomalous change of the equation of state. 

At present we have only an upper bound on tensor perturbations derived from COBE and WMAP, 
namely, H(4>{f = Ho)) < 2 x 10~ 5 Mp;. 9 ) Without taking a possible change in V[4>] during inflation 
into account, this corresponds to a constraint Qcwif ~ 1Hz, ao) < 8 x 10™ 16 . Thus WMAP data do 
not provide us direct information on the gravitational wave background. If we consider more specific 
inflationary models, however, we can predict a precise value of £Igw m the high-frequency region. 
Among many mechanisms of inflation, chaotic inflation proposed by Linde 17 ) is the most attractive 
for its naturalness of initial condition. WMAP, however, has ruled out one of them, namely a model 
with a quartic potential V[<p] = A(/? 4 /4. 6 ) The only remaining renormalizable model with a massive 
scalar potential, on the other hand, has been shown to be not only in good agreement with WMAP 
data 18 ) but also well motivated by particle physics and phenomenologically successful if we identify 
the inflaton with a sneutrino field. 19 ) This model predicts r = 0.16 on large scales, 18 ) which is 
within the reach of BICEP. 11 ) We can also calculate VLqw m the observable frequency range by 
laser interferometers as Qcw(f ~ lHz,ao) = 2.5 x 10~ 16 . Thus allowing a room of uncertainty in 
inflationary models, it is desirable to make a detector sensitive to 10~ 20 < Qgw < 10~ 15 . 

The sensitivity of a laser interferometer is often represented by the dimensionless gravitational 
wave amplitude h c (f). For a detection of the stochastic wave background, correlation analysis using 
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multiple detectors is very powerful, 1,2 ) and a detector with sensitivity h c (f) can probe Slew U P to 

where AT is the observational period, A/ is the effective bandwidth, and the coefficient X rep- 
resents the overlap of two detectors and is of order unity. 1 ' 2 ) The initial LIGO has a sensitivity 
corresponding to VL GW ~ 5 x 1CT 6 at / ~ 100Hz with AT = 4 months. 2 ) The advanced LIGO 
can detect Qgw ~ 5 x 1CT 11 for the same observational period. As is apparent from eq. (20), a 
lower frequency band is more advantageous for the detection of gravitational waves with a nearly 
scale-invariant spectrum, VLqw — constant. LISA, which is scheduled to be launched in 2011, has 
the best sensitivity at / ~ 10~ 3 Hz. 20 ) Although we cannot make correlation analysis with LISA due 
to the vanishing overlap of data streams whose noises do not correlate, we could attain the level of 
flew > 10~ 12 . 21 ) Unfortunately, however, the LISA band is severely contaminated by astrophysical 
sources and would hamper measurement of the primordial background with View < 10 -11 , because 
there are so many close white dwarf binaries that they cannot be resolved separately, which result 
in a confusion noise. 22 ) Their contamination is expected to disappear around / > 0.2Hz due to 
their large radii. The frequency band / ~ 1Hz is expected to be relatively clean, 21,23 ) and this 
band could be observed by DECIGO. 23 ) If its ultimately intended sensitivity is achieved, we could 
detect the stochastic background up to VLqw ~ 10~ 20 at ~ 0.2Hz. Thus DECIGO would serve as 
an ideal detector to study the stochastic gravitational wave background from inflation. 

The comoving wave corresponding to / ~ 1Hz today reentered the Hubble horizon around t ~ 
3 x 10~ 20 sec with the cosmic temperature T ~ 3 x 10 6 GeV in the standard cosmology. Thus, if 
significant entropy production takes place around t ~ 10~ 20 sec, we can detect its imprint on the 
spectrum of View using DECIGO. 

A typical mechanism of entropy production in this epoch is decay of massive particles, which we 
denote by %, with the decay rate T x ~ 10~ 7 — 10 _5 GeV. If this particle decays through gravitational 
interaction, its mass corresponds to m ~ (0.8 — 4) x 10 10 GeV. Evolution of the energy density of 
two components, the radiation (relativistic particles) and the particle x is solved with the following 
equations 

^^ = -±H Prad (t) + r xPx (t), *&& = -3H Px (t)-r xPx {t), (2i) 

where T x is the decay rate of the particle \ to radiation. These equations are solved as 

n f s ( a \ - 3 

PxW = Px\ a i 

-4 / „ \ -4 



- e - r *(^), a = a(t), (22) 



Prada) = /WK) (-) + T x (-) f ^l Px (a t )e- r ^-^dr, (23) 
where ai = a{ti) and U <C T" 1 is an initial time with p ra d{o.i) 3> p x (a«). 
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Figure 1 depicts log Qcwifi a o) m chaotic inflation model driven by a massive scalar field with 
various values of T x and the entropy increase factor F. The latter is defined by 

_ s(t e )a 3 (t e ) 4tt 2 3_4/w 

* - s(U)al ' S ~ 90 9 * 1 ~ 3 T ' 1 ' 

with s being the entropy density, and a* is the effective number of relativistic degrees of freedom. 
Here t e 3> T^ 1 is an arbitrary time in the radiation dominated regime after entropy production 
has terminated. In this figure, solid line represents the case with p x = from the beginning as a 
calibrator. The weak deviation from a flat straight line is due to the weak time dependence of the 
Hubble parameter during inflation as described by the second term in the right hand side of eq. 
(14). In the present model of a massive scalar inflaton, this curve is represented as &Gw(f,ao) = 
2.5 x 10 -16 [1 — 0.1 In (//1Hz)]. In this figure, dash-dotted line at the bottom depicts the case with 
T x = 10 _7 GeV with the entropy increase by a factor F = 10 3 . Three intermediate curves, on the 
other hand, all correspond to the case F = 10 but with different values of T x . 

As seen there, the ratio of the values of flew in the left and the right plateau depends only on 
the entropy increase factor (apart from the weak dependence of the solid line on /). The reason is 
simple. As the mode / reenters the Hubble radius, Qcwif,®) evolves as eq. (9) or 

M/> a) = ^^S Gwif, a m (f)), (25) 

with the initial condition given in eq.(8). Let us consider Q,aw{f, a) at a = a e = a(t e ) with 
Un(f) < t e . If U n (f) S> T" 1 the mode reenters the horizon after the entropy production has 
terminated, so the factor in the right hand side of eq. (25) is equal to unity and we find Qcw(f, o-e) = 
^Gw(f,ain(f)) in this case, which corresponds to the left plateau in figure 1. For the modes that 
reenter the horizon sufficiently early, we have ptot{flin)<Hn = Prad(o-in)<if n , so we find QcwifjO-e) = 
F~ A I^Q,Gw{fi ciin(f)), because ptot is also equal to p ra d at a = a e . Here we have assumed that g* 
does not change between U n (f) and t e . This region corresponds to the right plateau in the figure. 
Thus the ratio is given by F~ 4 / 3 apart from the weak dependence on V[4>(f)]. 

If the entropy production occurs at a later epoch with T x <C 10 -7 GeV, the entire frequency band 
observable with DECIGO lies in the right plateau region in the figure. We could still estimate the 
entropy increase factor F in such a case if we could measure the Hubble parameter during inflation 
with CMB polarization and extrapolate it to a higher frequency band. 

So far we have assumed that inflation-produced tensor perturbations are the only source of 
the stochastic gravitational wave background apart from possible contamination from Galactic and 
extra-Galactic stellar binaries. But, there are several proposed cosmological sources of backgrounds, 
such as cosmic strings 24 ) and field rearrangement associated with a global phase transition. 25 ) The 
former produces 

n GW (mz, a ) = 3 x 10' 8 ( 1Ql 6Q eV ) . (26) 
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where is the symmetry breaking scale of string formation. 26 ) This would exceed inflationary 
gravitational radiation if vl > 10 8 GeV. On the other hand, global phase transition induces 

n GW (f,a in (f))^Y^ (27) 

which should be compared with eq. (8). 25 -* Thus, in order to suppress the contribution related 
to the symmetry breaking scale, vg should be smaller than H((f>) during inflation. Even if these 
contributions surpass that of inflation, however, we can probe equation of state in the early universe 
with these gravitational waves because they also have intrinsically scale-invariant spectrum at the 
outset. (Note that eqs. (26) and (27) apply to the case phase transitions occur after inflation. If 
the symmetry breaking scales are so large that phase transitions take place in the early stage of 
inflation, they do not contribute to the stochastic gravitational wave background. Alternatively, if 
phase transitions occur in a late stage of inflation, 27 ) they might leave an interesting imprint on 
the stochastic background observable with laser interferometers.) 

In summary we have proposed to use the stochastic gravitational wave background from inflation 
to probe the equation of state of early universe and its evolution. Due to the nearly scale invariant 
nature of inflation, small deviations from a simple radiation dominated universe (w = 1/3) would 
be clearly imprinted in the spectrum of Slcwif-, a o)- We could obtain interesting information that 
can be hardly extracted by other methods. 

N.S. would like to thank A. Cooray for useful discussions and JSPS for a fellowship to re- 
search abroad. This work was supported in part by JSPS Grant-in-Aid of Scientific Research No. 
13640285(JY). 
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Fig. 1. Energy spectrum of the gravitational wave background, Q,Gw{f, do), in a chaotic inflationary model driven 
by a massive scalar field. Solid line represents the case with no entropy production. Dash-dotted line depicts the 
case with T x = 10 -7 GeV and the entropy increase factor F — 10 3 . Three intermediate curves are for F = 10 
with T x = 1CT 5 GeV (short-dashed line), T x = 1CT 6 GeV (dotted line), and T x = 1CT 7 GeV (long-dashed line), 
respectively. 



10 



